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Abstract. In this paper we introduce a working generalization of the theory 
of Grobner bases for algebras of partial difference polynomials with constant 
coefficients. One obtains symbolic (formal) computation for systems of lin- 
ear or non-linear partial difference equations arising, for instance, as discrete 
models or by the discretization of systems of differential equations. From an 
algebraic viewpoint, the algebras of partial difference polynomials are free ob- 
jects in the category of commutative algebras endowed with the action by 
endomorphisms of a monoid isomorphic to N r . Then, the investigation of 
Grobner bases in this context contributes also to the current research trend 
consisting in studying polynomial rings under the action of suitable symmetries 
that are compatible with effective methods. Since the algebras of difference 
polynomials are not Noetherian ones, we propose in this paper a theory for 
grading them that provides a Noetherian subalgebras filtration. This implies 
that the variants of the Buchberger's algorithm we developed for difference 
ideals terminate in the finitely generated graded case when truncated up to 
some degree. Moreover, even in the non-graded case, we provide criterions 
for certifying completeness of eventually finite Grobner bases when they are 
computed within sufficiently large bounded degrees. We generalize also the 
concepts of homogenization and saturation, and related algorithms, to the 
context of difference ideals. The feasibily of the proposed methods is shown 
by an implementation in Maple that is the first to provide computations for 
systems of non-linear partial difference equations. We make use of a test set 
based on the discretization of concrete systems of non-linear partial differential 
equations. 

1. Introduction 

An important idea at the intersection of many algebraic theories consists in 
studying algebraic structures under the action of operators of different nature, typ- 
ically automorphims and derivations. Classical roots of this idea can be found 
clearly in invariant and representation theory, as well as in the study of polynomial 
identities satisfied by associative algebras. Recently, topics like algebraic statistic 
[4] or entanglement theory [22] have given new impulse and applications to the 
research on such themes. Another fundamental source of inspiration is the theory 
of differential and difference algebras introduced in the pioneeristic work of Ritt 
[23l [24] and afterwards developed by Kolchin [16] , Cohn [6] , Levin [21] and many 
others. From the point of view of computational methods, starting from the algo- 
rithms proposed by Ritt himself, a considerable advancement can be recorded in 
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the differential case (see for instance [25]). Much less has been achieved for the 
algebras of difference polynomials where working algorithms can be found mainly 
in the linear case [T^. Nevertheless, the interest for such computations is relevant 
because of applications in the discretization of systems of differential equations like 
the automatic generation of finite difference schemes or the consistency analysis 
of finite difference approximations (9J [TTJ [19]. The present paper contributes to 
this research trend by concerning the development of effective methods for systems 
of linear or non-linear partial difference equations with constant coefficients. We 
provide also an implementation of such methods which is the first to allow com- 
putations in the non-linear case. Specifically, we generalize the theory of Grobner 
bases and related algorithms for ideals of the algebra of partial difference polyno- 
mials. We are able to do this in a general and systematic way, by defining classes 
of suitable monomial orderings, by extending the Buchberger's algorithm and the 
concept of grading to difference ideals, by defining truncated homogeneous compu- 
tations and even by introducing a suitable notion of homogenization for such ideals. 
First contributions to such theory can be found in [9] \17\ ITS]. In particular, owing 
to the notion of letterplace correspondence/embedding introduced in [17l [18], note 
that the Grobner bases computations for ideals of the free associative algebra are 
a subclass of the same computations for ideals of the algebra of ordinary difference 
polynomials. 

The algebras of partial difference polynomials are free algebras in the class of 
commutative algebras that are invariant under the action by endomorphisms of a 
monoid isomorphic to N r . Then, the study of Grobner bases for such algebras be- 
longs to the general investigation of computational methods for commutative rings 
or modules that have suitable symmetries. Moreover, from the viewpoint of appli- 
cations, the algebras of partial difference polynomials are fundamental structures 
in the formal theory of partial difference equations where a set of unknown mul- 
tivariate functions is assumed algebraically independent together with all partial 
shifts of them. To provide symbolic computation for systems of such equations is 
hence essential to introduce Grobner bases methods. Based on a suitable definition 
of monomial orderings that are compatible with shifts action and the description of 
large classes of them, the present paper introduces variants of the Buchberger's algo- 
rithm for partial difference ideals. These procedures take advantage of the monoid 
symmetry essentially by killing all S-polynomials in a orbit except for a minimal 
one. Note that the algebras of difference polynomials are not Noetherian since they 
are polynomial rings in an infinite number of variables and hence termination is not 
generally guaranteed for the proposed algorithms. With the aim of improving this 
situation, we define suitable gradings that are compatible with the monoid action 
and provide nitrations of the algebra of partial difference polynomials with finitely 
generated subalgebras. We obtain therefore the termination for finitely generated 
graded difference ideals when computations are performed within some bounded 
degree. For non-graded ideals but for monomial orderings compatible with such 
gradings, we prove also criterions able to certify that a Grobner basis computation 
performed over a suitable finite set of variables that is within a sufficiently large de- 
gree, is a complete one. Finally, the paper generalizes the notion of homogenization 
and saturation to difference ideals with respect to the given gradings and provides 
the algorithms to perfom this ideal operations. As a byproduct, one obtains an 
alternative algorithm to compute Grobner bases of non-graded difference ideals via 
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homogeneous computations. By means of an implementation in Maple, all these 
methods are finally experimented on difference ideals obtained by the discretization 
of systems of non-linear differential equations. 

2. Algebras of partial difference polynomials 

Fix K any field and let E be a monoid (semigroup with identity) that we denote 
multiplicatively. Let A be a commutative if-algebra and denote End^-(A) the 
monoid of .fT-algebra endomorphisms of A. We call A a ^-invariant algebra or 
briefly a Ti-algebra if there is a monoid homomorphism p : E — > Endx(A). In this 
case, we denote a ■ x = p(a)(x), for all a £ E and x £ A. Let A,B be E-algebras 
and tp : A — > B be a if-algebra homomorphism. We say that <p is a "E-algebra 
homomorphism if tp(a ■ x) ~ a ■ tp(x), for all a £ E and x £ A. Let A be a E-algebra 
and let I C A be an ideal. We call / a Yi-invariant ideal or simply a Yi-ideal if 
S • J C J. Clearly, all kernels of S-algebra homomorphisms are E-ideals. 

Definition 2.1. Let A be a 'S-algebra and let X d A be a subset. We say that A 
is E-generated by X if A is generated by E • X as K-algebra. In other words, A 
coincides with the smallest E-subalgebra of A containing X. In the same way, one 
defines Yi-generation for the E-ideals. 

In the category of E-invariant algebras one can define free objects. In fact, let 
X be a set and denote x(a) each element (x, a) of the product set X(E) = X x E. 
Define P = K[X(E)] the polynomial algebra in the commuting variables x(a). For 
any element a £ E consider the if-algebra endomorphism a : P — >• P such that 
x(t) i — y x(crr), for all x(t) £ ^(E). Then, one has a faithful monoid representation 
p : E — > End^-(P) such that p(a) = a ad hence P is a E-algebra. Note that if E is 
a left-cancellative monoid then all maps p(o~) are injective. 

Proposition 2.2. Let A be a T,-algebra and let f : X — s> A be any map. Then, 
there is a unique Yi-algebra homomorphism tp : P A such that <p{x(l)) = f{x), 
for all x £ X . 

Proof. It is sufficient to define ip(x(a)) = a ■ f(x), for all x £ X and a £ E. In fact, 
one has tp(r ■ x(o~)) = (p(x(ra)) = ra ■ f(x) = r • (a ■ f{x)) = r • ip(x(o~)), for any 

r e E. □ 

Definition 2.3. We call P — if[X(E)] the free E-algebra generated by X . In fact, 
P is E- generated by the subset ^(1). 

In other words, the algebra P is an essential tool in the theory of E-algebras 
because any such algebra A that is E-generated by a set X can be obtained as a 
quotient E-algebra P/I, where / is a E-ideal of P. For instance, from the viewpoint 
of computational methods, if one develops them for P then such methods can be 
extended to any quotient P/I as it is done in the classical theory of Grdbner bases 
for affine algebras. Note also that if E is defined as the monoid Inc(N) = {/ : N — >• 
N | / strictly increasing}, or some power of this, one obtains an environment for 
computations in algebraic statistic [1]. 

We want now to go in the direction of developing fundamental structures for sym- 
bolic (formal) computation on systems of partial difference equations with constant 
coefficients. From now on, we assume that X = {xq, x%, . . .} is a finite or countable 
set and E is a free commutative monoid generated by a finite set, say {eri, . . . , o>}. 
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Then, we consider the free E-algebra P = K[X(T,)}. Note that (E, •) is a cancella- 
tive monoid isomorphic to (N r ,+) and the monomorphisms p(a) : P — > P have 
infinite order for all er 7^ 1. For any Xi(o~) G ^(S), we call i and a respectively the 
index and the weight of the variable Xi(o~). If we put X(o~) = {xi(o~) \ Xi G X} and 
Si(S) - [Xiia) I (t £ E} one has clearly P = ® aeS K[X{a)} = ® x . eX K[ Xi (E)], 
where all subalgebras K \X(a)] are isomorphic to K[X] and all subalgebras if[xi(E)] 
toK[E]. 

Definition 2.4. The free ^-algebra P = K[X(E)} (E = (a 1} . . .,oy}) is called the 
algebra of partial difference polynomials with constant coefficients. 

The motivation for such name is in the formal theory of partial difference equa- 
tions [6l [21]. In this theory, in fact, the indeterminates Xi(l) are by definition 
algebraically independent unknown functions . . . , t r ) in the variables tj and 

the maps p{(Jk) are the shift operators . . . , t r ) t-> Ui(t\, . . . , tk + h, . . . ,t r ) 

where ft is a parameter (mesh step). If a — Y[i a T z then the indeterminates 
Xi((j) = <y ■ Xi(l) correspond to the (algebraically independent) shifted functions 
Ui(t\ + a±h, . . . , t r + a r h) — a ■ Uj(ti, . . . , t r ). Then, a E-ideal / C P is also called a 
partial difference ideal and a E-basis of I corresponds to a system of partial differ- 
ence equations in the unknown functions itj(ii, . . . , t r ). One uses the term ordinary 
difference when r = 1. Note that the algebras of difference polynomials are not 
Noctherian rings since they are polynomial rings in an infinite number of variables. 
One has therefore that difference ideals have bases or E-bases which are generally 
infinite. 

In the next sections we generalize the Grobner basis theory to the free E-algebra 
P = K [X(E)] of partial difference polynomials. Clearly, one reobtains the classical 
theory when E = {1} (r = 0) that is P — K\X]. The starting point is to define 
monomial orderings of P which are compatible with the action of the monoid E. 

3. Monomial E-orderings 

Denote by M = Mon(P) the set of all monomials of P. Note that even if the set 
X(E) is infinite (in fact countable), one can endow P by monomial orderings. This 
is an important consequence of the Higman's Lemma |15] which can be stated in 
the following way (see for instance pQ, Corollary 2.3 and remarks at beginning of 
page 5175). 

Proposition 3.1. Let -< be a total ordering on M such that 

(i) \ <m for all m G M; 

(ii) -< is compatible with multiplication on M , that is if m -< n then tm -< tn, 
for any to, n, t G M. 

Then -< is also a well-ordering of M that is a monomial ordering of P if and only 
if the restriction of -< to the variables set X(H) is a well-ordering. 

Clearly, it is easy to construct well-orderings for the set X(H) which is in bijective 
correspondence to N r+1 . Note that the monoid E stabilizes the monomials set M 
since it stabilizes X(E). We introduce then the following notion. 

Definition 3.2. Let -< be a monomial ordering of P. We call -< a (monomial) 
E-ordering of P if -< is compatible with the H-action on M , that is m -< n implies 
that a ■ to -< a ■ n for all m,n G M and o G E. 
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A straightforward consequence of this definition is the following result. 

Proposition 3.3. Let -< be a monomial E-ordering of P. Then m < a ■ m for all 

m G M and a G E. 

Proof. By contradiction, assume that there are ra, a such that m >- cr ■ m. Then, 
a ■ m >- a 2 ■ m and by induction one obtains the infinite descending chain m >- 
a ■ m >- a 2 ■ m >- . . . which contradicts that -< is a well-ordering. □ 

The orderings on the variables set X(E) that can be extended to monomial 
E-orderings are as follows. 

Definition 3.4. Let -< be a well-ordering o/X(E). We call -< a (variable) E- 
ranking of P if -< is compatible with the T,-action on X(E), that is u -< v implies 
that a ■ u -< cr • v /or aZ/ u, v G X(E) and ueS. 

As for Proposition 13. 3[ we have that if -< is a E-ranking then u ^ a ■ u for all 
u e AT(S) and a £ S. Moreover, if X is a finite set then condition u ^ cr • u for 
all u, cr together with E-compatibility implies that -< is a well-ordering by applying 
Dickson's Lemma (or Higman's Lemma) to E which is isomorphic to N r . However, 
note that in this paper the set X may be also countable. 

Owing to the decompositions X(E) — \J ae - s X(a) = \J x . eX x i(^) °^ tne var i" 
ables set of the ring P, we can define E-rankings of P in a natural way. Denote by 
Q the monoid if -algebra defined by the free commutative monoid E = (cr 1; . . . ,o>). 
In other words, Q — K[o~±, . . . ,o~ r ] is the polynomial algebra in the commutative 
variables crj . From now on, we assume that E is endowed with a monomial ordering 
< of Q. By abuse, we call < a monomial ordering o/E. 

Definition 3.5. Fix < a monomial ordering o/E. For all Xi(a),Xj(r) G X(E), 
we define: 

(i) Xi (cr) -< Xj (r) if and only if a < r or a — r and i < j . In other words, 
X(a) -< X(t) when a < r. 

(ii) Xi{o~) -<' Xj(r) if and only if i < j or i = j and a < r . In other words, 
cCi(E) -<,' Xj(E) when i < j. 

Clearly -< and -< are 6oi/i "E-rankings of P that we call respectively weight and 
index E-ranking defined by a monomial ordering of E. 

For all x t G X and a G E denote P(ct) = if [X(a)], M(a) = Mon(F(cr)) and 
P{xi) = K[xi(E)], M(xi) — Mon(P(xi)). Owing to the decompositions P = 
^gj. P(c) = (£) x eX P{xi), one has that a monomial m G M can be uniquely writ- 
ten as m = m(di) ■ ■ ■ m(dk) — m(xi 1 ) ■ ■ ■m(xi l ), where m(5 p ) G M(S p ),m(xi ) G 
M{xi p ) and c?i > . . . > 6k, i± > . . . > ij. By means of such presentations we can 
define block monomial orderings of P extending weight and index ranking. Recall 
that p : E — 5- End/<(P) is the faithful monoid representation defined by the action 
of E over P. For any a G E one has that the map p(cr) defines an isomorphism 
between the monoids M(l),M(cr) and hence between the algebras P(l),P(cr). In 
other words, we have M(cr) = a ■ M(l), P(c) = cr • P(l). 

Definition 3.6. Fix -< a monomial ordering of the subalgebra P(l) C P and extend 
it to all subalgebras P(cr) (cr G E) by the isomorphisms p(o~). In other words, we 
put a ■ m -< cr • n if and only if m -< n, for any m,n G M(l). Then, for all 
m,n G M, m — m(Si) ■ ■ ■ m(Sk), n = n(8\) ■ ■ ■ n(8k) with 5± > . . . > 6k we define 
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to -< w n if and only ifm(5j) — n(6j) if j < i and m(Si) -< n(Si) for some 1 < i < k. 
Clearly, the restriction of -< w to the variables of P is just the weight Ti-ranking. 

Proposition 3.7. The ordering -< w is a T,-ordering of P. 

Proof. Note that if to = m(Si) ■ ■ ■ m(Sk) G M with m(<Ji) G M(<7j) and 6% > 
. . . > Sk then a ■ m — m(a8\) ■ ■ • m(a5k), where m(aSi) = a ■ m(Si) G M{a5i) and 
aS\ > . . . > aSk since < is a monomial ordering of E. Assume to -< w n that is 
m(Sj) = n(5j) for j < i and m(5i) -< n(#i). Clearly m(aSj) = n(o~5j) for j < i and 
one has m(#,) -< n((5i) if and only if to(1) -< n(l) if and only if m(aSi) -< n(cr5i). 
Then, we conclude that er • m ~< w a ■ n. □ 

Note that we have also a monoid faithful representation : N — >• Endi<-(P) 
such that the endomorphism </>(i) is defined as Xj(o~) H > Xi+j(cr) for any i,j > 
and cr € E. Clearly 0(z) induces isomorphism between the monoids M(xq), M(xi) 
and the algebras P(xq), P(xi). The algebra P(xo) can be easily endowed with a E- 
ordering. For instance, since P(xq) = ® CTgS K[xo(a)} one can define a lexicographic 
ordering as in Definition 13.61 

Definition 3.8. Fix -< a monomial T.-ordering of the subalgebra P(xq) C P and 
extend it to all subalgebras P{xi) (Xi G X) by the isomorphisms cf>(i). For any 
to, n G M, to = 771(2;^ ) ■ • ■ m(a;, fc ), n — n{xi t ) ■ ■ ■ n(xi h ) with i\ >...>%}. we put 
to -<j n if and only if m(xi ) = n{xi ) if q < p and m(xi p ) -< n(xt p ) for some 
1 < p < k. Note that the restriction of -<j to t/ie variables of P is the index 
Yi-ranking. 

Proposition 3.9. The ordering -<i is a T.-ordering of P. 

Proof. Note that if m = m(xi 1 ) ■ ■ ■m(xi h ) G M with m(xi ) G M(xi p ) and ii > 
. . . > i k then a ■ m = m'ix^) ■ ■ ■m'(x ik ) where m'{xi p ) = a ■ m{x ip ) G M{x lp ). 
Suppose to -<j n that is m{xi q ) = n(xi ) if 9 < p and m(xi p ) -< n(xi p ). We have 
clearly that m'(xi ) = n'(xi ). Moreover, since -< is a E-ordering of P(xq) and 
therefore of P(xi p ), one has also m'(xi p ) -< n'(xi p ) that is cr • to -<j cr • n. □ 

We call the above monomial E-orderings of P respectively weight E- 

ordering defined by a monomial ordering of P{\) and index Yi-ordering of P defined 
by a monomial T,-ordering of P(xq). Clearly, both these orderings depend also on 
a monomial ordering of E. Note that index E-orderings are suitable for generation 
of finite difference schemes for partial differential equations [101 E] • The weight E- 
orderings are instead compatible with the gradings of the E-algebra P we introduce 
in Section 5. For this reason they are suitable for obtaining complete Grobner bases 
from partial computations. 

4. Grobner E-bases 

From now on, we consider P endowed with a monomial E-ordering -<. Let 
/ = J2i c i m i ^ P with TOj G M, Ci G K, a ^ 0. We denote as usual lm(/) = — 
max^{m,}, lc(/) = c k and lt(/) = lc(/)lm(/). If G C P we put lm(G) = {lm(/) | 
/ G G, f ^ 0} and we denote as LM(G) the ideal of P generated by lm(G). 

Proposition 4.1. Let G C P. Then lm(E • G) = E • lm(G). In particular, if I is 
a H-ideal of P then LM(7) is also H-ideal. 
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Proof. Since P is endowed with a E-ordering, one has that lm(cr • /) = a ■ lm(/) for 
any JeF,//0 and a € £. Then, £ • lm(7) = lm(£ • I) C lm(7) and therefore 
LM(7) = (lm(J)) is a S-ideal. □ 

Definition 4.2. Let I C P be a H-ideal and G C I. We call G a Grobner E-basis 
of I iflm(G) is a Yi-basis o/LM(J). In other words, £ ■ G is a Grobner basis of I 
as P -ideal. 

Since the monoid £ is assumed isomorphic to N r that is S-ideals are partial 
difference ideals, we may say that Grobner E-bases are partial difference Grobner 
bases [9]. Another possible name is H-equivariant Grobner bases [4]. Simplicity and 
generality lead us to the previous definition that already appeared in |f 8j . 

Let f,g S P, f,g ^ and put lt(/) = cm, \t(g) = dn with m,n £ M and c, d G K . 
If I = lcm(m, n) we define as usual the S-polynomial spoly(/, g) — {I / cm) f ~ {I / dn)g . 
Clearly spoly(/, g) = -spoly(g, /) and spoly(/, /) = 0. 

Proposition 4.3. For all f,g G P,f,g ^ and for any a G £ one has a ■ 

spoly(/, g) = spoly(cr ■ f,a ■ g). 

Proof. Since £ acts on the variables set A(£) by injective maps, it is sufficient to 
note that a ■ lcm(m, n) — lcm(<7 • m,a ■ n) for all m,n G M and a G £. □ 

The following definition is a standard tool in Grobner bases theory. 

Definition 4.4. Let f G P, / ^ and G G P . // / = with f l G P, g l G G 

and lm(/) y lm(fi)lm(gi) for all i, we say that f has a Grobner representation 
respect to G. 

Note that if / = Yli f*9i 1S a Grobner representation then a- f = J2i( a ' fi)( a ' 9i) 
is also a Grobner representation, for any cr G £. In fact, since -< is a E-ordering 
of P one has that lm(/) >^ \m(fi)\m(gi) implies that lm(<7 • /) = a ■ lm(/) >z 
(a ■ hn(/i))(<j • lm(gi)) = lm(cr-/i)lm(cr-f/i) for all i. A celebrated result from Bruno 
Buchberger [$: is the following. 

Proposition 4.5 (Buchberger's criterion). Let G be a basis of the ideal I G P. 
Then, G is a Grobner basis of I if and only if for all f,g G G, f, g ^ the S- 
polynomial spoly(/, g) has a Grobner representation with respect to G. 

Usually the above result, see for instance [8], is stated when P is a polynomial 
algebra with a finite number of variables and G is a finite set. In fact, such as- 
sumptions are not needed since Noetherianity is not used in the proof, but only 
the existence of a monomial ordering for P. See also the comprehensive Bergman's 
paper [2 t where the "Diamond Lemma" is proved without any restriction on the 
finiteness of the variables set. We want now to prove a generalization of the Buch- 
berger's criterion for Grobner E-bases of P. For this purpose it is useful to introduce 
the following notations. 

Definition 4.6. Let a = JJi^i' S T = IL °f e s - We denote gcd(a,r) = Ut^T 
where ji — min(aj, fti) , for any i. 

Proposition 4.7 (E-criterion). Let G be a Y^-basis of a H-ideal I G P. Then, 
G is a Grobner T,-basis of I if and only if for all f,g G G,f,g ^ and for any 
a, t G E such that gcd(cr, r) = I , the S-polynomial spoly(er • /, r • g) has a Grobner 
representation with respect to E • G. 
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Proof. We prove that E-G is a Grobner basis of I and we make use of the Proposition 
14.51 Then, consider any pair of elements a-f, r-g 6 E-G where /, <? £ G, /, <? ^ and 
a, t £ E. Put (5 = gcd(<7, r) and hence cr = <5(t', r = St' with cr', r' £ £, gcd(cr', r') = 
1. By Proposition l4~3l we have spolvfcr- f,r-g) — <5-spoly((r'-/, r'-g). By hypothesis, 
assume that spoly(er' • /, r' • g) — h — J2 U fv{v • g^), with v £ £, /„ £ P, g„ £ G, is 
a Grobner representation with respect to E • G. Since -< is a E-ordering of P . we 
conclude that we have also the Grobner representation spoly(<r ■ /, r ■ g) — 5 ■ h — 
E«V ■ f v ){5v ■ g v ). ' q 

A standard procedure in the Buchberger's algorithm is the following. 



Algorithm 4.1 Reduce 
Input: G C P and / £ P. 

Output: he P such that f -he (G) and h = or lm(/i) £ LM(G). 

while /i ^ and lm(/i) £ LM(G) do 

choose g e G, g ^ such that lm(g) divides lm(/i); 
ft := h - (\t(h)/lt(g))g; 

end while; 

return h. 



Note that the termination of Reduce is provided since -< is a monomial ordering 
of P. In particular, even if G is an infinite set, there are only a finite number of 
elements g £ G,g ^ such that lm(g) divides \m(h) and hence hn(g) ^ lm(ft). 
It is well-known that if Reduce(/, G) = then / has a Grobner representation 
with respect to G. Moreover, if Reduce(/, G) = h ^ then clearly one has 
Reduce(/, G U {h}) — 0. Therefore, from Proposition 14.71 it follows immediately 
the correctness of the following algorithm. 



Algorithm 4.2 SigmaGBasis 

Input: H, a E-basis of a E-ideal I C P. 
Output: G, a Grobner E-basis of /. 
G := H; 

B:={(f,g)\f,g€G}; 
while B ^ do 

choose £ B; 

B:=B\{(f,g)}; 

for all a, t £ E such that gcd(er, t) = 1 do 
h := REDUCE(spoly(cr •/, r •<?),£ • G); 
if /i ^ then 

B:=BU{( Jl /i) 1 (/» l l»)| ff 6G}i 
G:=GU{fe}; 
end if; 
end for; 
end while; 
return G. 
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Clearly, all well-known criteria (product criterion, chain criterion, etc) can be 
applied to SigmaGBasis to shorten the number of S-polynomials to be considered. 
In fact, one can understand this algorithm as the usual Buchberger's procedure 
applied to the basis E • H , where an additional criterion to avoid "useless pairs" is 
given by Proposition [+7j Owing to Non-Noetherianity of the ring P, note that the 
termination of SigmaGBasis is not provided in general and this is, in fact, one 
of the main problems in differential/difference algebra. Nevertheless, in the next 
section we introduce some suitable grading for the algebra P which provides that 
a truncated version of the algorithm SigmaGBasis with homogeneous input stops 
in a finite number of steps. Some variant of the algorithm SigmaGBasis appeared 
in [9] and before in [17j [18] for the ordinary difference case. 

5. An illustrative example 

In this section we apply the algorithm SigmaGBasis to a simple example in 
order to provide a concrete computation with it. Let X = {x, y}, E = (<ri, 02) and 
consider the algebra of partial difference polynomials P = K[X(E)] . To simplify the 
notation, we identify the monoid (E, •) with (N 2 , +) by means of the isomorphism 
a i a 2 l— * Then, we denote the variables of P as x(i, j),y(i, j), for all i,j > 0. 

We consider now the E-ideal (difference ideal) I C P that is E-generated by the 
difference polynomials 

51 =i/(l, 1)2/(1, 0)-2x(0,l) 2 , 

52 =2/(2,0)+; Sf(0, 0)1(1,0). 

In other words, this E-basis (difference basis) encodes a system of non-linear differ- 
ence equations with constant coefficients in two unknown bivariate functions. By 
symbolic (formal) computations, we want to substitute this system with a comple- 
tion of it, namely a Grobner E-basis. We may want to do this for the purposes of 
checking membership of other equations to the E-ideal, elimination of unknowns, 
etc. The main problem is that such basis may be infinite, but it is not the case for 
this example. We fix then the degrevlex ordering on E {a\ > 02) and the lex mono- 
mial ordering on K[x,y] {x >~ y) that is on P(0, 0) = K[x(0, 0), y(0, 0)]. A weight 
E-ordering (Definition I3-6|) is hence defined on P and we use it for computing a 
Grobner E-basis of I. Such basis consists of the elements 51 , 52 together with the 
difference polynomials 

53=y(l,2M0,l) 2 -y(l,0)x(0,2) 2 , 

54 = 2.t(1, l) 2 - i(0, 0)a;(l, 0)x(0, l)x(l, 1). 

Note that in all these elements the first monomial is the leading one with respect 
to the given monomial ordering of P. Let us see how the algorithm SigmaGBasis 
is able to obtain such Grobner E-basis. Since the E-ideal I is E-generated by 
G = {51, 52} then / is generated as an ideal of P by E • G that are the polynomials 

(m) ' .91 = V(i + 1, 3 + + " MhJ + I) 2 ; 
(hj) -92 =y(i + 2,j) + x(i,j)x(i + l,j), 

for all i,j > 0. By applying the product criterion, we have to consider three kinds 
of S-polynomials 

spoly((i,j + 1) ■ 51, ■ gi) = (i, j) ■ spoly((0, 1) -51,51), 
spoly((i + 1, j) • 51, -52) = ■ spoly((l,0) -51,52), 
spoly((i + 1, j + 1) • 51, (i, j + 2) • 52) = (i, j + 1) • spoly((l, 0) • 51, (0, 1) • 52). 
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The E-criterion implies therefore that one has to reduce with respect to the basis 
E • G just the S-polynomials 

si = spoly((0,l) -51,51), 
s 2 = spoly((l,0) -51,52), 
S3 -spoly((l,0)-5i,(0,l)-5 2 ). 

The reduction of the S-polynomial si leads to the new element 53 and the current 
E-basis of / is now G — {51,52,53}- The additional S-polynomials that survive to 
product and E-criterion are 

s 4 = spoly((0,l) -51,53), 

s 5 = spoly((0,2) -51,53), 

s 6 = spoly((0,2)-5 2 ,(l,0)-5 3 ). 
We have that S4 — > and S2 —> 54 with respect to E ■ G. The E-basis is then 
G = {51,52,53,54} and one has a new S-polynomial 

s 7 = spoly((l,0) -53,54)- 

Finally, we have that all S-polynomials S3, S5, sq, S7 reduce to zero with respect to 
S • G and hence G is a Grobner E-basis of /. Note that G is in fact a minimal such 
basis and also that in this simple example there is no use of the chain criterion that 
can be always applied together with the other criteria. 

6. Gradings of P compatible with E- action 

We want now to introduce some gradings of the algebra P = K[X(E)] which 
are compatible with E-action and formation of least common multiples in M — 
Mon(P). As before, we fix a monomial order < of E. We start extending the 
structure (E,max, ■) in the following way. 

Definition 6.1. Let be an element disjoint by E and put E = E U {0}. Then, 
we define a commutative idempotent monoid (E, +) with identity that extends the 
monoid (E, max) (with identity 1) by imposing that + er = a, for any a £ E. More- 
over, we define a commutative monoid (E, ■) with identity 1 extending the monoid 
(E, •) by putting 0-cr = 0, for all a 6 E. Since multiplication clearly distributes over 
addition, one has that (E, +, •) is a commutative idempotent semiring, also known 
as commutative dioid |13j . 

Note that the faithful monoid representation p : E — > End^ (P) can be extended 
to E where p(0) : P — > P is the algebra endomorphism such that Xi(a) M> 0, for all 

XiWei(s). 

Definition 6.2. Let w : M — > E be the unique mapping such that 

(i) w(l) = 0; 

(ii) w(mn) — w(m) + w(n), for any m,n G M ; 

(iii) w(xi(a)) = a , for all i > and a € E. 

Note that (i),(ii) state that w is a monoid homomorphism from the free commutative 
monoid (M, ■) to (E, +). We call w the weight function of P. 

More explicitely, if m = Xi 1 (5i) ai ■ ■ ■Xi k (5k) ak is any monomial of P different 

from 1 then w(m) = Si H h 6k = max<((5i, . . . , 5k) ■ We denote M a = {m £ M \ 

w(m) = a} and define P a C P the subspace spanned by M a , for any a £ E. Because 
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w : (M, •) — > (E, +) is a monoid homomorphism one has that P = © CTe f; Pa is a 
grading of the algebra P over the commutative monoid (£, +). If / G P a we say 
that / is a w -homogeneous element and we put w(/) = a. Recall that for any a G E 
we denoted P(o") = -ftr[X(<j)] which is a subalgebra of P = (E)] isomorphic to 
If we put P(0) = P =K then one has that P^ = r < ff P CT = (g) T < (J P(r) 
is a subalgebra of P. In particular, we have that P^ — Po © Pi = P(0) ® P(l) = 
P(l) is isomorphic to the polynomial algebra 

Definition 6.3. A monomial order < o/E is said sequential if {t G E | t < er} is 

a finite set, for all a G E. 

It is important to note that if X is a finite set and < is a sequential ordering 
of E then the sequence {P^ a G E} is a filtration of P consisting of Noether- 
ian subalgebras. For such reason, from now on we assume E be endowed with a 
sequential monomial ordering. 

Proposition 6.4. The weight function satisfies the following properties: 

(i) w(er ■ m) = erw(m), /or any cr G E and m G M; 

(ii) w(lcm(m, n)) = w(mn) — w(m) + w(n), for all m,n G M. Then, m \ n 
implies that w(m) < w(n). 

Proof. If m = 1 then w(<7 • m) = w(m) = = crw(m). If otherwise m = 
aJi! {o~i) ai ■■■x ik (Sk) ak with (5i > . . . > 4 then a ■ m = x h (a5i) ai ■ ■ ■ x ik (oS k ) ak 
where <7<5i > . . . > aS k since < is a monomial ordering of E. We conclude that 
w(er • m) = a5\ = aw(m). To prove (ii) it is sufficient to note that the weight of a 
monomial does not depend on the exponents of the variables occuring in it. □ 

Note that the property (i) implies that the map w is a homomorphism with 
respect to the action of E on M and E. In other words, one has that aP T C P aT 
for any a G E,r G E. Moreover, the property (ii) means that w is also a monoid 
homomorphism from (M, lem) to (E,+). 

Definition 6.5. Let I be an ideal of P. We call I a w-graded ideal if I — I a 
with l a = I n Pj. In this case, if I is also a Y,-ideal then a ■ I T C I aT for all 
a G E,t G E. 

Owing to the w-grading of P, one can show that a truncated version of the algo- 
rithm SigmaGBasis admits termination. If f,g G P, / ^ g are w-homogeneous el- 
ements then the S-polynomial h = spoly(/, g) is clearly w-homogeneous too. More- 
over, by property (ii) of Proposition 16.41 we have that w(h) — w(/) + w(g) and 
hence if w(f),w(g) < S then also w(h) < 5, for some S G E. By means of this 
remark, one obtains immediately the following result. 

Proposition 6.6 (Truncated termination over the weight). Let I C P be a w- 

graded Y,-ideal and fix S G E. Assume I has a w-homogeneous basis H such that 
Hg = {/ G H | w(/) < 5} is a finite set. Then, there is a w-homogeneous Grobner 
H-basis G of I such that G is also a finite set. In other words, if we consider 
for the algorithm SigmaGBasis a selection strategy of the S -polynomials based on 
their weights ordered by <, we obtain that the S -truncated version of SigmaGBasis 
stops in a finite number of steps. 

Proof. First of all, note that the algorithm SigmaGBasis computes essentially a 
subset G of a Grobner basis E • G obtained by applying the Buchberger's algorithm 
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to the basis E • H of /. Moreover, by Proposition l6.4l the elements of E • H and hence 
of E • G are all w- homogeneous. Denote H' s = {<r • / | a € E, / € H, crw(/) < 5}. 
Since < is a sequential monomial order of E and H$ is a finite set one has that 
H' s is also a finite set. We consider therefore Xg the finite set of variables of P 
occurring in the elements of H' s and define P^s) = K[Xg] C P. In fact, the 6- 
truncated algorithm SigmaGBasis computes a subset of a Grobner basis of the 
ideal Irg^ C P(s) generated by H' s . By Noetherianity of the finitely generated 
polynomial ring P^ we clearly obtain termination. □ 

Clearly the above result provides algorithmic solution to the ideal membership 
for finitely generated w-graded E-ideals. Note that if r = that is E = {1} then 
the algorithm SigmaGBasis coicides with classical Buchberger's algorithm and 
Proposition 16.61 states that if / is a finitely generated ideal of P = Po ® Pi = 
K[xq, x\, . . .] then I has also a finite Grobner basis. According with the above 
proof, this is a consequence of the fact that the Buchberger's algorithm runs over 
the finite number of variables occuring in the generators of /. 

Another useful grading of P can be introduced in the following way. Consider 
the set N = NU{- oo} endowed with the binary operations max and +. Then 
(N, max, +) is clearly a commutative idempotent semiring (or commutative dioid 
or max-plus algebra). Define deg : E — > N the mapping such that deg(O) = — oo 
and deg(er) = a*, for any a = ]X erf' . Clearly deg is a monoid homomorphism 
from (£,•) to (N,+). 

Definition 6.7. Let ord : M — > N be the unique mapping such that 

(i) ord(l) = — oo; 

(ii) ord(mn) = max(ord(m), ord(n)), for any m,n € M; 
(hi) ovd(xi(a)) = deg(er), for all i > and a G S. 

Clearly (i),(ii) state that ord is a monoid homomorphism from (M, •) to (N, max). 
We call ord the order function of P. 

For any monomial m = a:^ (ch) Ql • ■ • X{ k (Sk) ak different from 1 we have that 
ord(m) = max(deg((5i), . . . , deg(<$fc)). Clearly, the order function defines a grading 
P = (J) de ^ Pd 01 the algebra P over the commutative monoid (N, max). Define 
P (d) = ®i< d p i = ®dc g ( CT )< ( i P ( CT ) which is a subalgebra of P. Then, if X is 
a finite set we have that the sequence is a filtration of P with 

Noetherian subalgebras where P^ = -P_oo © Pq is isomorphic to 

Definition 6.8. A monomial order < of E is said compatible with deg when 
deg(cr) < deg(r) implies that a <t, for any a, r G E. 

If < is compatible with deg, note that < is a sequential ordering of E and 
ord(m) = deg(w(m)), for all m G M. Finally, one has that the weight and order 
functions clearly coincide when r = 1. 

Proposition 6.9. The order function verifies the following: 

(i) ord(cr • m) = deg(er) + ord(m), for any a G E and m G M ; 

(ii) ord(lcm(m, n)) = ord(mn) = max(ord(m), ord(n)), for allm,n G AI . There- 
fore, if m | n then ord(m) < ord(n). 

Proof. For m = 1 one has ord(er ■ m) = ord(m) = — oo = deg(cr) + ord(m). If oth- 
erwise m = Xi ± (Si) ai ■ ■ ■ Xi k (Sk) ak then a-m = Xi t (aSi) ai ■ ■ ■ xi k (a5k) ak and hence 
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ord(cr • m) = max(deg(cr(Si), . . . , deg(aS k )) = deg(a) + max(deg(6i), . . . , deg(S k )) = 
deg(cr) + ord(rn). Property (ii) follows immediately as in Proposition ^. 41 □ 

Definition 6.10. Let I be an ideal of P. We call I a ord-graded ideal if I = J2i h 
with Ii ~ I n Pi. If I is also a T,-ideal then a ■ Ii C Idcg(a)+i f or an V °~ €E E and 
i G N. 

Consider now /, g G P, f ^ g ord-homogeneous elements. The S-polynomial h = 
spoly(/, g) is clearly ord-homogeneous and ord(/i) = max(ord(/), ord(<?)). Then 
ord(/), ord(g) < d implies that ord(ft.) < d, for some d G N and one proves the 
following result as for Proposition 16. 61 

Proposition 6.11 (Truncated termination over the order). Let I C P be a ord- 
graded Yi-ideal and fix d G N. Assume I has a ord-homogeneous basis of H such 
that Hd — {/ G H | ord(/) < d} is a finite set. Then, there is a ord-homogeneous 
Grobner "E-basis G G I such that Gd is also a finite set. In other words, if we 
consider for SigmaGBasis a selection strategy of the S- polynomials based on their 
orders, we have that the d-truncated version of SigmaGBasis terminates in a finite 
number of steps. 

By means of weight and order functions one has criterions, also in the non-graded 
case, that provide that a Grobner E-basis is the eventually finite complete one even 
if it has been computed within some bounded weight or order for the algebra P 
that is over a finite number of variables. This is of course important because actual 
computations can be only performed in such a way. As before, we fix a sequential 
monomial ordering < on E. 

Definition 6.12. Let -< be a monomial E- ordering of P. We call -< compatible 
with the weight function ifw(m) < w(n) implies that m -< n, for all m,n s M. In 
a similar way, one defines when -< is compatible with the order function. 

Proposition 6.13. Let ~< w be a weight ^-ordering as in Definition \3.b\ Then -< w 
is compatible with the weight function. In particular, if the monomial order < of E 
is compatible with deg then -< w is also compatible with the order function. 

Proof. Let m — m{5\) ■■■ m(Sk), n — n(S±) ■ ■ ■ n{8k) two monomials of P with 
m(5i),n(Si) € M(Si) and d\ > . . . > 5k- Assume m -< w n that is m(5j) — n(Sj) if 
j < i and m(Si) -< n(5i) for some 1 < i < k. If i > 1 or m(Si) =/= 1 then clearly 
w(m) = w(n) = 8±. Otherwise, we conclude w(m) < 5i = w(n). Moreover, if < is 
compatible with deg then ord(m) = deg(w(m)) < deg(w(n)) = ord(n) implies that 
w(t7j) < w(n) and hence m < w n. □ 

For any S G E,d G N define now E 5 = {a G E | a < 6} and E d = {a G E | 
deg(a) < d}. 

Proposition 6.14 (Finite E-criterion). Assume the Yi-ordering of P be compatible 
with the weight function. Let G C P be a finite set and denote I the Ti-ideal 
generated by G. Moreover, define 8 — max<{w(lm(g)) | g G G}. Then, G is a 
Grobner T,-basis of I if and only if for all f,g£G and for any a, t G E such that 
gcd(er, t) = 1 and gcd(<rdm(/), r-lm^)) ^ 1, the S-polynomial spoly(cr- f,r-g) has a 
Grobner representation with respect to the finite set E^2 • G. In the same way, if the 
Ti-ordering of P is compatible with the order function and d — max{ord(lm(g)) | 
g G G}, then G is a Grobner Yi-basis of I when the above S-polynomials have a 
Grobner representation with respect to Y,2d • G. 
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Proof. Let spoly(cr • /, r • g) = h = ^ y f v (v ■ g v ) be a Grobner representation with 
respect to E • G that is \m(h) y lm(f v )(u ■ lm(g v )) for all v. We want to bound 
the elements v G E with respect to the ordering <. Put m = lm(/),n = Im(g) 
and hence lm(<7 ■ /) = a ■ m, lm(<7 • g) — a ■ n. By product criterion, we can 
assume that u = gcd(er • m, r ■ n) ^ 1. Then, there is a variable ^(cra) = Xi(rf3) 
that divides it where ^(a) divides m and hence a < w(m) < 5 and £,(/3) divides 
n and therefore j3 < w(n) < (5. Then era = r/3 and one has that a \ ft, r | a 
because gcd(er, r) = 1. We conclude that a, r < S and if v — lcm(<7 • m,r • m) 
then w(w) = max(crw(m), rw(n)) < <5 2 . Clearly v >- lm(/i) ^ v ■ lm((;^) and 
hence S 2 > w(v) > vw{\ai{g l/ )) > v. In a similar way, one argues for the order 
function. □ 

The above criterion implies that with respect to S-orderings compatible with 
weight or order functions one has an algorithm able to compute a finite Grobner 
S-basis, whenever this exists, in a finite number of steps. In practice, this results in 
an adaptative procedure that keeps the bound S 2 (2d) for the elements of £ to be 
applied to the generators, constantly updated with respect to the maximal weight 
5 (order d) of the leading monomials of the current generators. In particular, note 
that the bound is S (d) when forming S-polynomials and S 2 (2d) for reductions. 

7. Homogenizing with respect to order function 

The purpose of this section is to analyze (de)homogenization processes in the 
context of S-ideals. Such methods are generally developed to have structures and 
computations that arc homogeneous with respect to some grading, even if the input 
data are not such. Besides to the theoretical advantages as the concept of projective 
closure, these techniques usually imply computational advantages (see for instance 
[3]). Note that for E-ideals one can develop these methods for both w-grading and 
ord-grading, but we will concentrate on the second one. Univariate homogenizations 
are usually more efficient than multivariate ones because leading monomials are 
preserved by the homogenization process. 

Let t be a new variable disjoint by X. Define X = X U {£},X(£) = X x 
S,P = i4TLY(£)] and finally M = Mon(P). Consider the algebra endomorphism 
tp : P P such that Xi(a) i-> xi(a) and t(a) H> 1, for all i,a. Clearly <p 2 = ip and 
P = <p(P). Moreover, one has that ip is a E-algebra endomorphism. Then <p defines 
a bijective correspondence between all E-ideals of P and E-ideals of P containing 
ker<^ = (t(l)-l) s . 

Definition 7.1. Denote by N = iV or d the largest ord-graded H-ideal contained in 
ker ip that is the ideal generated by all oid-homogeneous elements f € P such that 
Vif) = o. 

Proposition 7.2. The ideal N C P is generated by the elements 

(i) t(a) — t(r) for all a, r G E, a ^ r, deg(er) = deg(r); 

(ii) t(a)t(r) — t(a),x(a)t(r) — x(a) for any a, t G E, deg(cr) > deg(r). 

Proof. Let / G P be a ord-homogeneous element such that <p(f) — 0. Since the 
polynomials of type (i),(ii) clearly belongs to N, we have to prove that / is congruent 
to modulo (i),(ii). Assume first that all variables of / belong to £(S) = {t(cr) \ a G 
E}. Recall that if m = t(Si) ai ■ ■ ■ t(b~k) ak is any monomial of / then d = ord(/) = 
max(deg(5i), . . . , deg((5/c)). Therefore, one has that / is congruent modulo (ii) to 
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/' = X^ C ^( T where n £ E,deg(ri) = d and a £ K,J2i c i = 0- B y applying 
identity (i) it follows that / is congruent to Ci)t(a) = for some fixed a such 
that deg(cr) = d. 

Consider now the general case when the variables of / belong to X(E). Fix 
a £ E such that deg(er) = d. Modulo the identities (i),(ii), one has that / is 
congruent to a polynomial /' whose monomials are either of type m £ M such that 
ord(rn) = d or of type t{a)n where n £ M, ord(n) < d. We show that in fact /' = 0. 
Denote /' = t(a)g — h where g, h are polynomials in P, h is ord-homogeneous and 
ord(/i) = d. Since = ip(f') = g — h one has that /' = (t(a) — l)g. If we 
assume j / then the monomials n of g are such that ord(n) = d which is a 
contradiction. □ 

We want now to define a bijective correspondence between all E-ideals of P and 
some class of ord-graded E-ideals of P containing N. 

Definition 7.3. Let I be any H-ideal of P. We define I* C P the largest ord- 
graded Y,-ideal contained in the preimage ip^ 1 ^) that is I* is the ideal generated 
by all ord-homogeneous elements in <p~^(I). Clearly N = 0* C I*. We call I* the 
ord-homogenization of the Ti-ideal I. 

Definition 7.4. Let f £ P, f ^ and denote f — fd the decomposition of 
f in its ord-homogeneous components. We define topord(/) = d' = max{<i}. If 
f £ K that is d' = —oo we put f* = f. Otherwise, we denote f* — t(a)f where 
a £ E such that deg(cr) = d! . We call topord(/) the top order of f and f* its 
ord-homogenization. 

Proposition 7.5. Let I be a E-ideal of P. Then I* = (f* \ f £ I,f ^ 0) + N. 

Proof. Denote J = (/* | / e I) + N. Clearly J C I*. Let g £ I* be a ord- 
homogeneous element and define / = <p(g) £ I. If / = then g £ N C J. 
Otherwise, denote d = topord(/) and d' — ord(g). Since clearly d 1 > d one has 
that g is congruent modulo N to h = t(a)f, where a £ E such that deg(cr) = d 1 . 
Hence, if d! = d then h is congruent exactly to /*. Otherwise, the polynomial h is 
congruent to t(a)f* . In both cases, we conclude that g is congruent modulo N C J 
to an element of J and therefore g £ J. □ 

If I C P is a E-ideal one has clearly that = L. Moreover, if J C P is a 

ord-graded E-ideal containing N then in general J C (p(J)* ■ 

Definition 7.6. Let N C J C P be a oid-graded H-ideal. Define J' = f(J)* — 
I / € J, f ^ N, f ord-homogeneous) + N. Then J C J' C P is a ord-graded 
Ti-ideal that we call the saturation of J. 

Definition 7.7. Let J C P be a ord-graded 12 -ideal containing N . We say that J is 
saturated if J coincides with its saturation ip(J)* that is for any ord-homogeneous 
element f £ J, f ^ N one has that <£>(/)* £ J. If I is a T,-ideal of P then its 
ord-homogenization I* is clearly a saturated ideal. 

Therefore, a bijective correspondence is given between all E-ideals of P and the 
saturated ord-graded E-ideals of P containing N. 

We want now to analyze the behaviour of Grobner E-bases under homogenization 
and dehomogenization. Note that the arguments of Proposition 17.21 implies clearly 
that the polynomials (i),(ii) are in fact a Grobner basis of the ideal N with respect 
to any monomial ordering of P. For this reason we introduce the following notion. 
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Definition 7.8. A monomial m G M is said normal modulo TV if m G M or 
m = t(a)n with n G M, a G £ such that d — deg(er) > ord(n). Moreover, we 
assume that t(a) = min^{<(r) | deg(r) = d}. A polynomial f G P is in normal 
form modulo TV if all its monomials are normal modulo TV. 

Definition 7.9. Let -< be a £- ordering of P compatible with the order function. 
We call -< a ord-homogenization E-ordering if t(a)m -< n /or aZZ m, n G M, tr G S 
sitc/i </iai deg(er) = ord(n) > ord(m). 

It is easy to define one of the above orderings. Fix for instance the lex or 
degrevlex monomial order on the polynomial ring P(l) — K[xq(1), xi(l), . . . , t(l)] 
where xo(l) > Xi(l) > ... > t(l). Moreover, fix a monomial ordering on £ which 
is compatible with deg and define the weight £-ordering ~< w of P as in Definition 
13. 61 Clearly -< w is a ord-homogenization S-ordering. 

From now on, we assume P be endowed with a ord-homogenization £-ordering. 

Proposition 7.10. Letp,q G M be two normal monomials modulo N such that 
ord(p) = ovd(q). Then p < q implies that (p(p) -< f{q)- 

Proof. By definition, the monomials p, q are of type m G M or t(o~)m with deg(cr) > 
ord(m). Since -< is a ord-homogenization order, when comparing two of such mono- 
mials of the same order one has only the following cases: m -< n, t(o~)m -< t(o~)n 
or t(a)m -< n. Then, we have to prove tp(p) = m -< n = only when 

i(er)m -< n. This follows immediately from -< is compatible with the order function 
and ord(m) < ord(n) = deg(cr). □ 

From now on, for any / G P. f ^ we denote by /* the normal form of t(a)f 
modulo TV where a G E, deg(er) = topord(/). 

Proposition 7.11. Let /ePJ/O be a oid-homogeneous polynomial in normal 
form modulo N. Then lm(<^(/)) = (^(lm(/)). Moreover, we have that lm(/*) = 
lm(/) forallfePJ^O. 

Proof. The first part of the statement follows immediately from Proposition 17.101 
Moreover, if a G E,deg(cr) = topord(/) then t(a) cannot appear in the leading 
monomial of /* and hence lm(/*) = lm(/). □ 

Definition 7.12. Let N C J C P be a Yl-ideal. Moreover, let G C J be a subset 
of polynomials in normal form modulo N. We say that G is a Grobner E-basis of 
J modulo N if G U TV is a Grobner Ti-basis of J. 

Proposition 7.13. Let TV C J C P be a ord-graded 'E-ideal. Lf G is a ord- 
homogeneous Grobner Yi-basis of J modulo TV then (p(G) is a Grobner Yi-basis of 
<p(J). 

Proof. Since G is a Grobner E-basis of J modulo TV we have that for any ord- 
homogeneous polynomial / G J, / 7^ in normal form modulo TV there is an element 
g G G and a G E such that a ■ lm(g) | hn(/). Then, by applying the E-algebra 
endomorphism <p one obtains that a - \m(ip(g)) | lm(<p(/)) that is <p(G) is a Grobner 
E-basis of ip(J). □ 

Proposition 7.14. Let L C P be a Y^-ideal and let G be a Grobner Y.-basis of L. 
Then G* = {g* | g G G} is a oid-homogeneous Grobner basis of I* modulo TV. 
Moreover, one has that lm(G*) = lm(G). 
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Proof. Let /' G I* be a ord-homogencous clement in normal form modulo N and 
put / = ip(f) e I. Then, either /' = /* or /' = t(a)f* with ord(/') = deg(cr) > 
ord(/*) = topord(J). Since G is a Grobner E-basis of / there is g G G, t G E 
such that r • lm(g) | lm(/). By Proposition 17. 1 II one has that lm(/) = lm(/*) and 
lm(<7) = lm(g*). Therefore, r • \m(g*) divides lm(/*) and this monomial clearly 
divides lm(/'). □ 

By the above propositions we obtain immediately what follows. 

Corollary 7.15. Let N C J C P be a ord- graded Y.-ideal and denote J 1 = <p(J)* 
its saturation. Moreover, let G be a ord-homogeneous Grobner Ti-basis of J modulo 
N . Then G' = <p(G)* = {(f{g)* \ g G G} is a ord-homogeneous Grobner Yi-basis of 
J' modulo N . Moreover, we have lm(G') = \m(ip(G)). 

Let / C P be any S-ideal. The previous results suggest an alternative method 
to calculate a Grobner E-basis of I which is based only on ord-homogeneous com- 
putations. Assume H is any E-basis of I and denote as before H* = {/* | / S H}. 
Clearly J = + N is a ord-homogeneous E-ideal of P containing N such that 

<p(J) = I. Assume now we compute G a ord-homogeneous Grobner E-basis of J 
modulo N. Then, ip(G) is a Grobner E-basis of /. Note that by using a ord-based 
selection strategy for the S-polynomials, the Grobner E-basis G can be obtained 
order by order automatically minimal that is a ■ hn(/) not divides lm(g) for all 
f,g £ G, f ^ g and a G E. This is clearly a computational advantage, but since 
generally lm(G) ^ lm(y>(G)) one has that ip(G) may be not minimal. In the worst 
case, the ideal J may have an infinite and hence uncomputable minimal Grobner 
E-bases but / just a finite one. This is clearly not the case when one considers a 
saturated ideal J' = I* since we have lm(G') = \m(ip(G')) when G' is a minimal 
Grobner E-basis of J'. Note that this nice property depends on the fact that we 
deal with a univariate homogenization. A drawback is that if one computes the 
saturation J' by means of the ideal J according to Corollary 17. 151 one has again to 
compute a Grobner E-basis of J. Then, a better approach consists in computing 
"on the fly" the Grobner E-basis of J' starting from the generating set {/* | / G H}. 
In other words, any time that a new generator g of the ord-homogeneous Grobner 
E-basis arises from the reduction of an S-polynomial, we saturate g that is we sub- 
stitute this polynomial with ip(g)*. In formal terms, the algorithm one obtains is 
the following one. 
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Algorithm 7.1 SigmaGBasis2 

Input: H, a E-basis of a E-ideal I C P. 

Output: ip(G), a Grobner E-basis of I such that lm(G) = lm(ip(G)). 
G:=H*; 

B:={(f,g)\f,geG}; 

while B 7^ do- 
choose (/, g) e B; 
B:=B\{(f,g)}; 

for all ct, r e S such that gcd(er, t) = 1 do 
/i := REDUCE(spoly(tr • /, r • g), E • G U AT); 
if /i ^ then 

ft = <p{h)* 

B:=BU{(g,h),(h,h)\geG}; 
G :=GU{h}; 
end if; 
end for; 
end while; 
return f{G). 



Proposition 7.16. The algorithm SigmaGBasis.2 is correct. 

Proof. Note that at each step we are inside an ideal J such that ip( J) = I that is 
whose saturation is J' = I* . Moreover, for any ord-homogeneous element h G P 
which is in normal form modulo N one has that h' — (p(h)* divides h. This implies 
that if an S-polynomial is reduced to zero by adding h to the basis G, the same 
holds if we substitute h with h' . In case of termination, owing to the set G is a ord- 
homogeneous Grobner E-basis of J modulo N whose elements are all saturated, by 
Corollary 1 7 . 1 5 1 we may conclude that J — J' and hence <p(G) is a Grobner E-basis 
of / such that lm(G) = lm(>(G)). □ 

About termination or just termination up to some order d, this is not provided 
in general for the above algorithm. The reason is that even if all computations are 
ord-homogeneous, because of the saturation h = <p(h)* that may decrease the order 
we can't be sure at some suitable step that we will not get additional elements of 
order < d in the steps that follow. 

8. An illustrative example (continued) 

We apply now the algorithm SigmaGBasis2 to the same E-ideal that has been 
considered in Section 5 for illustrating SigmaGBasis. Recall that such ideal is 
I = {91,92) C P = K[X(E)] where X = {x,y},E = N 2 and 

<7i =j/(l,l)j/(l ) 0)-2a:(0,l) 2 
Sa = 1,(2,0) + *((], 0)ar(l,0). 

Let now X = {x,y,t} and define the polynomial algebra P — ATLY(E)] with vari- 
ables x(i,j),y(i,j),t(i,j), for all i,j > 0. We consider the E-algebra endomorphism 
<p : P —> P such that x(i,j) h-> x(i,j),y(i,j) n> y(i,j) and t(i,j) h-> 1. In Propo- 
sition EH we proved that the largest ord-graded E-ideal contained in kercp is the 
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ideal N C P generated by the polynomials 

t(i, j) — t(k, I) where (i, j) ^ (k, I), i + j = k + l; 

t(i,j)t(k, I) — t{i, j), x(i,j)t(k, I) — x(i, j) where i + j > k + I. 

Moreover, we define a ord-homogenization E-ordering of P (Definition 17. 9p as the 
weight E-ordering given by the degrevlex ordering of E and the lex monomial or- 
dering of K[x, y,t] (x >- y >- t). Note that such E-ordering extends the one defined 
for P C P in Section 5. Recall that a Grobner E-basis of / with respect to such 
ordering is given by the elements 51,52 together with 

ff3 =y(l,2M0,l) 2 -y(l,0):r(0,2) 2 , 

54 = 2x(l, If - x(0, 0)o;(l, 0)a;(0, l)o;(l, 1). 

We introduce then ord-homogenizations of the polynomials g\ , gi that are 

ff t = y(l,l)y(l,0)-2*(0,2)x(0,l) 2 , 
g* 2 =y(2,0)+t(0,2)x(0,0)x(l,0). 

We start applying SigmaGBasis2 to G — {gl, g%} by considering the S-polynomials 

si = spoly((0,l) -g$,g$), 
s 2 = spoly((l,0) -Sfijffa), 
.s 3 = spoly((l,0)- 5l *,(0,l)- 52 *). 

By reducing s± with respect to E • G U N one obtains exactly the polynomial 

gl = y(l, 2)x(0, l) 2 - t(0, 3)y(l, 0)x(0, 2) 2 . 

Clearly ^3 is already a saturated element. Then G = {g* 1 g\ , g^} and we form the 
new S-polynomials 

s 4 = spoly((0,l) -gt,gt), 
s 5 = spoly((0,2)- 5l *,^), 
S6 = spoly((0,2)- 5 *,(l,0)- ff3 *). 
With respect to E • G U N, one has the reductions S4 — > and s 2 h where 

/i = 2t(Q, 3)x(l, l) 2 - i(0, 3)x(0, 0)x(l, 0)x(0, l)x(l, 1) = t(0, 3)g|. 

By saturating this element as ip(h)* — gl, we update G = {gl , g 2 , g%, g%} and 
another S-polynomial is defined as 

s 7 = spoly((l,0) -53,54). 

All remaining S-polynomials S3, S5, sq, S7 reduce to zero with respect to S-GUiV 
and we conclude that <p(G) — {51,52,53,54} is a Grobner E-basis of / C P. 

9. Testing and timings 

In this section we present a set of tests for the algorithms SigmaGBasis and 
SigmaGBasis2 which is based on an experimental implementation of them in the 
language of Maple. This is actually the first implementation of algorithms for com- 
puting Grobner bases of systems of linear or non-linear partial difference equations. 
Note that for the linear case one has the packages LDA (Linear Difference Algebra) 
[12"] and Ore_algebra[shift_algebra] in the Maple distribution. The main idea that 
lead us when coding the proposed algorithms is that they can be considered variants 
of the classical Buchberger's algorithm where some amount of computations can be 
avoided by means of the symmetry defined by the monoid E. In fact, as explained in 
the previous sections, a "basic" approach to calculate a Grobner E-basis of a E-ideal 
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/ generated by a E-basis H consists in applying the Buchberger's algorithm to the 
basis E • H . One obtains therefore a Grobner basis G of / from which a Grobner 
E-basis G C G can be extracted such that E • lm(G) = lm(G'). Clearly, chain 
and coprime criterions can be used in the usual way in the procedure. Then, the 
algorithm SigmaGBasis can be understood as the variant that prescribes the appli- 
cation also of the E-criterion (Proposition ^. 7|) to the S-polynomials spoly(<7-/, r -g) 
and to add the set of all shifts E • h to the current basis when a new element h 
arises from the reduction of an S-polynomial. Then, the Grobner E-basis of / is 
simply the union of the initial basis H with the new elements h. Recall that the 
procedure is correct only if one uses a monomial E-ordcring. Clearly, from the set 
E is infinite it follows that actual computations can be only performed with a finite 
subset of E that is over a finite set of variables of P = K[X(T,)]. Typically, one 
fixes a bound d for the degree of the elements of E that is for the order of the 
variables Xi(a). Owing to the finite E-criterion f Proposition l6. 14]) . a basis obtained 
with a monomial ordering compatible with the order function is certified to be a 
complete Grobner E-basis if the order bound is at least the double of the maximum 
top order of its elements. 

In addition to the basic procedure for the computation of Grobner E-bases and 
the algorithm SigmaGBasis, for the experiments we consider also a variant of the 
latter method where the E-criterion is suppressed but one continues to shift the 
reduced form of the S-polynomials. This procedure is tested to the aim of under- 
standing the contribution of any of the implemented strategies. Finally, we propose 
an implementation of the algorithm SigmaGBasis2 based on the saturation of a 
E-ideal with respect to the grading defined by the order function. In practice, once 
one has homogenized the initial generators, the saturation (p(h)* is performed be- 
fore the application of shifting, for each new element h obtained by the reduction 
of an S-polynomial. In output one returns the dehomogenization of the computed 
basis. Note that this procedure is correct only if one uses a E-ordering which is 
compatible with the order function and if the polynomials are kept in normal form 
modulo N during the computations. 

The monomial E-orderings of P that we consider for the tests are defined in 
the following way. One has initially to fix a monomial ordering for E and we 
choose degrevlex in order to provide compatibility with the degree. Then, one fixes 
a monomial ordering, for us lex, over the subring P(l) = _?£T[a"(1)] or P(xq) — 
K[xq(Y>)} that is extended as a block ordering to the polynomial ring P = K[X(E)] 
according to the choice of a variables ranking based on weight or index respectively. 
We distinguish these two cases in the examples by the letters "w" and "i". The 
integer that comes before these letters refers to the fixed order bound. Note that 
the algorithm SigmaGBasis2 is compatible only with rankings of type weight. 

For the basic variant of the Grobner E-bases algorithm, one can clearly use any 
implementation of the Buchberger's algorithm as, for instance, the one contained 
in the package Groebner of Maple. We have preferred instead to develop ourselves 
all different variants in order to have the same implementation and hence the same 
efficiency, for the fundamental subroutines of the algorithms. In this way, for the 
basic version we have been also able to access to important parameters of the com- 
putation as the total number of S-polynomial reductions. This number is for us the 
sum of the actual S-polynomials with the initial generators that are interreduced. 
Note that our implementation of the Buchberger's algorithm is in fact generally 
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comparable with the built-in ones of Maple. For instance, the test falkow-6w-basic 
takes 9 hours, but using Groebner[Basis] it takes 11 hours with method=buchberger 
and 8.5 hours with method=maplef4. Other parameters that are considered for the 
experiments are the number of input and output generators. Note that for the basic 
algorithm we count generators and not ^-generators. Finally, the parameter "mi- 
nout" refers to the number of elements of a minimal Grobner S-basis. All examples 
have been computed with Maple 12 running on a server with a four core Intel Xeon 
at 3.16GHz and 64 GB RAM. The timings are given in hour-minute-second format. 
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We give now some details about the examples we have used. All the examples are 
non-linear systems of ordinary or partial difference equations with constant coeffi- 
cients which are of interest in literature. For instance, the tests falkow are obtained 
by the discretization of the Falkowich-Karman equation which is a non-linear two- 
dimensional differential equation describing transonic flow in gas dynamics. The 
discretization we used are equations (41) in | TTj . Then, the navier examples are 
based on equations ei,e2,e3,e4 of the system (13) in the paper [TU] that are a fi- 
nite difference scheme corresponding to the discretization (9) of the Navier-Stokes 
equations for two-dimensional viscous incompressible fluid flows. The tests heat are 
the discretization of the one-dimensional heat equation as described in the equa- 
tions (10) and (11) of [20]. Finally, eq26 and eq27 are the equations (2.6) and (2.7) 
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at page 24 of [14] which are examples of ordinary difference equations that have 
periodic solutions. 

By analyzing the experiments, it is sufficiently clear that the strategy imple- 
mented in SigmaGBasis is the safest one and hence on the average, the most 
efficient one. In fact, by decreasing the number of S-polynomials this strategy 
avoids the dramatical effects of involved reductions as for the tests falkow-6w and 
navier-8w. For simpler examples the four strategies appear essentially equivalent. 
The algorithms SigmaGBasis and SigmaGBasis2 lead to practically identical 
computations but the latter method suffers of some overhead which is probably 
due to our still experimental implementation. For instance, even if the normal 
form modulo the ideal N is described in the Definition I7.8[ in our implementation 
we obtain it by computations that is by adding a Grobner basis of N to the input 
basis for SigmaGBasis2. 

The proposed algorithms usually provide only partial informations about the 
structure of Grobner E-bases since they are in general infinite. Nevertheless, it is 
interesting to note that by means of the finite E-criterion we have been able to 
certify that the examples falkow, navier and heat have finite bases with respect to 
the weight ranking. In particular, the elements of the Grobner basis of falkow have 
maximum top order equal to 4 and hence they are certified in order 8 in about 4 
minutes. The example navier has max top order equal to 6 and its certification is 
obtained in order 12 in less than one hour. Finally, the example heat has max top 
order 2 and it gets certification in order 4 in seconds. 



10. Conclusions and future directions 

This paper shows that one can not only generalize in a systematic way the 
Grobner bases theory and related algorithms to the algebras of partial difference 
polynomials but also make these methods really work by introducing suitable grad- 
ings for such algebras. In fact, weight and order functions provide a Noetherian 
subalgebras filtration that implies termination and completeness certification for 
actual computations that are performed within some bounded degree that is over 
a finite number of variables. We have then developed the first implementation of a 
variant of the Buchberger's algorithm for systems of linear or non-linear partial dif- 
ference equations. Even if such implementation is just experimental, the approach 
corresponding to the algorithm SigmaGBasis is strong enough to let it able to 
work with discretizations of real world systems of non-linear differential equations. 

In this paper we consider difference equations with constant coefficients and 
hence a next step along this line of research is to extend the proposed methods 
to systems of difference equations with non-constant coefficients that is to assume 
that E acts on the base field K in a non-trivial way. Moreover, since the algebras 
of partial difference polynomials are free objects in the category of commutative 
algebras endowed with the action of a monoid E isomorphic to N r , a natural future 
direction consists in extending the ideas introduced here to other types of monoid 
symmetry over commutative algebras as the ones used, for instance, in algebraic 
statistic [4]. Starting from Grobner bases, classical directions are the computation of 
the Hilbert series and free resolutions that one may generalize to partial difference 
ideals or other types of invariant ideals. Finally, we aim to have the proposed 
algorithms implemented in the kernel of computer algebra systems in order to tackle 
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involved problems related with the discretization of systems of partial differential 
equations [9l[TfllT2]. 
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